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The large impact parameter (bt) behavior of the solution is discussed as well as the 
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Glass Condensate) kinematic domain are analyzed. We obtain the asymptotic behavior 
for the physical observables and found the unitarity bounds for them. 
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1 Introduction 



High density QCD PQI21IH] which is dealing with the parton systems where the gluon occupation 
numbers are large, has entered a new phase of it's development: a direct comparison with the 
experimental data. A considerable success jU El El El El CHI HI] has been reached in description of 
new precise data on deep inelastic scattering [12] as well as in understanding of general features 
of hadron production in ion- ion collision J3] • The intensive theoretical f3J E3 E3 work lead to 
effective theory in the high parton density region with non-linear evolution equation for dipole- 
dipole amplitude that governs the dynamics in this region jTTj. The practical application 
of the high density QCD approach to a comparison with the experimental data is based on 
analytic [THI EH 120] and numerical El EE EH 122] solution to the non- linear equation. In spite 
of understanding of the main qualitative and, partly, quantative properties of the non-linear 
dynamics, this equation has not been solved at arbitrary values of the impact parameters (bt) 
and a number of "ad hoc" ansatzs were used for &i-dependence. 

It turns out that ^-dependence is a challenging problem since the non-perturbative correc- 
tions are important in the large b t kinematic region. There are two different ideas on the market 
how to include the non-perturbative large bt behavior into high density QCD dynamics. The 
first one [2H1 123 EH] claims that the non-perturbative corrections could be taken into account 
only in initial condition for the non-linear evolution equation, while the other idea demands 
the change of the kernel in the equation [26J. 

The objective of this paper is to find the semi-classical solution to the non-linear evolution 
equation with a special attention to the large b t behavior of the solution. The semi-classical 
approach has several advantages [T| ITi] II 5j. Firstly, it gives simple, transparent and analytic 
solution to the non-linear equation. Secondly, this solution has a good theoretical accuracy in 
the most interesting kinematic region: at low x and small sizes of interacting dipoles. Thirdly, 
this approach leads to a natural definition of the new saturation scale which is the principle 
dimensional parameter that governs dynamics at high density QCD region P3I21E]- It appears 
in the semi-classical approach as a critical line which divides the whole set of trajectories in 
two parts: inside and outside of the saturation region. In other words, this critical line is a 
separation line between two phase: color gluon condensate (CGC) and the gluon liquid and the 
saturation scale is the order parameter for this phase transition. 

The shortcoming of the semi-classical approach is the fact that the experimental data are 
not in the kinematic region where this approach has a good theoretical accuracy. However, we 
believe, that we can develop the reasonable method for the numerical solution of the non-linear 
equation only after studying the property of the semi-classical approach. 

In the next two sections we discuss the semi-classical approach for the dipole sizes smaller 
than the saturation scale (r t < 1/Q s (x) where Q s (x) is the saturation momentum). We review 
the semi-classical solution which has been investigated in this region (see Refs. PI El E3) 
concentrating mostly on large 6 r behavior of this solution. 

Section 4 is devoted to the semi-classical approach inside of the saturation region (r t > 
1/Q s (x)). We suggest a semi-classical solution which describes the behavior of the dipole- 
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dipole amplitude deeply inside of this domain and discuss the matching of this solution with 
the semi-classical solution at small dipole sizes. It turns out that this matching occurs on the 
special trajectory of the non-linear equation which defines the saturation scale for dipole-dipole 
rescattering. 

In section 5 we discuss the unitarity bound as well as energy behavior of the gluon structure 
function and the 7* — 7* total cross section. 

Section 6 gives the estimates of the accuracy of our approach. We calculate the first enhanced 
diagram as well as corrections to the non-linear equation that we neglected in our solution. 

In the last section we summarize our results and discuss the possible practical applications. 



2 General Approach 

2.1 The non- linear evolution equation 

The non-linear equation that governs the dynamics in the high parton density QCD domain 
can be written in the following form given by Balitsky and Kovchegov jTT] : 

dN(x 01 ,y;b t ) 2C F a s ( ATI , x , 2C F a s ( . 

In \ N(x 01 ,y;b t ) + x (2.1) 



' f X 2 / 1 1 

/ g? 2 x 2 2 01 2 2 iV(x 02 , y; b t - -X12) - iV(x 2, y\ b t - -xi 2 ) iV(x 12 , y\ b t 

. Jp XnoXio \ 1 Z 




The meaning of Eq. (j2.1j) is very simple and can be seen in Fig. ^ and Fig. 121 It describes 
the process of dipole interaction as two stages process. The first stage is decay of the initial 
dipole with size 201 into two dipole with sizes xyi and 202 with probability 



1 2 X Q\ 



|*(^oi -> ^02 + 3:12)1 - 2 2 • 

X 02 X \2 

In the second stage two produced dipoles interact with the target. The non-linear part of 
Eq. (|2.1jl takes into account the Glauber corrections for such an interaction ( see Fig. Q). The 
first term in Eq. ()2.1|) stands for possibility for the initial dipole to interact with the target 
without decaying into two dipoles. 

The equation is written in the coordinate representation which has a certain advantage since 
N(xoi,y]b t ) = Ima(x i,y; b t ) where a is the dipole amplitude. Directly from the unitarity 
constraint follows that N < 1 giving the natural asymptotic behavior for N, namely iV tends 
to unity at high energies. On the other hand the non-linear term in Eq. (|2.1jl contains the 
integration over sizes of produced dipoles. It means that we cannot conclude that the non- 
linear term is small even when the sizes of the initial dipoles are very small. 
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Non - Linear colour dipole evolution equation 




N(x 01 ) N(x 12 ) N(x 02 ) N( X 12) ■ N(x 02 ) 



Figure 1: The non-linear evolution equation. 



Emission term in momentum representation 



\ 



dJVjxoi, b t , y) 
dy 



where 



1C f ol s /x?.\ 



2 C F a s 



I"' 



x 2 





- _ZV(x 2, b t , y} 


X 2 X 2 
-*-02 ""-12 







JV(x 01 ,b t ;y) = lm a dipo i e (x 01 , b t ; y) 




Gluon Reggeization 
in momentum 
representation 



Probability for 
dipole -> 2 dipoles 
decay 



Figure 2: The linear (BFKL) equation. 
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It turns out to be useful to consider Eq. (|2.1|) in a mixed representation, fixing b t but 
introducing transverse momenta as conjugated variables to dipole sizes. The relation between 
these two representations are given by the following equations (see also Fig. |2J 1 \ 



oo 



N(x,y;b t ) = x 2 k dk J (kx) N(k,y;b t ) ; (2.2) 

JO 

poo Arr 

N(k,y;b t ) = / —J (kx)N(x,y;b t ) . (2.3) 
Jo x 

In momentum representation the nonlinear term in Eq. (|2.1|) . being the convolution in the 
coordinate representation, transforms to the product of two amplitudes, (N), at the same value 
of the initial transverse momentum (k). Therefore, the large value of the initial transverse 
momentum guarantees the smallness of the non-linear corrections. One can actually see that 
the non- linear term in Eq. (|2.1|) is a convolution only at large value of the impact parameter b t . 
Indeed, let us rewrite the non-linear term in Eq. (|2.1j) in the momentum representation going 
also to momentum transfer instead of bt- 

In this case 



(2vr 



N(k,y;Q) = I ^\e~^ N(k,y;b t ) . (2.4) 



The non-linear term reduces to the form 



d x 2 ^ — 2-iV(xo2, y; b t - -xi 2 )iV(xi2, y; b t - 77X02) 

X 02 X 12 ^ * 



J d 2 Q d 2 Q' 8(k -k' + ^Q- )-Q') N(k, y; Q) N(k', y; Q') e®'^ . (2.5) 

At large values of bt 3> 1/k and 1/k' both Q and Q' are small and of the order of 1/bt- 
Neglecting terms \Q and \Q' in 5-function in Eq. ()2.5|) we can reduce this equation to the 
product of N(k, y; b t ) ■ N(k, y; b t ). 

However, one can see that there is a region of integration over Q , namely, Q + Q' ~ l/b t 
k or/and k', while Q — Q' > k or/and k'. In the coordinate representation (see Eq. (|2.1jl ) this 
region of integration corresponds to X\i ~ X02 ~ 2 6 t ^> x i- I n this region the non-linear 
term can be reduced to 

/ d 2 ^N 2 (2b t ,y;Ab t ) (2.6) 

with Ab t = b t — |^02- In the region of large 6 t r 2 > x i = r\ (see Fig. Elfor all notations) 
the dipole amplitude in Eq. (J2.6|) describes the interaction of the dipole with the size which is 
much larger than the size of the lower dipole in Fig. |3J Such interaction cannot be described 
by Eq. (j2.1j) . Indeed, the non-linear equation is proven only in the case when r\ = xqi r 2 
since only in such conditions we can restrict ourselves by consideration of the "fan" diagrams 
(see Fig. EJ-a) and only in the case of the running QCD coupling [1. Eq. (j2.1j) is also valid for 

1 The extra factor x 2 in Eq. Ij2.2|l makes N dimensionless. 
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deep inelastic processes with nuclei [T7| but we will not consider this process here since in this 
case the dipole amplitude does not depend on b t for b t < Ra where Ra is the nucleus radius 
2 . Therefore, we should consider this region of integration as a correction to the non-linear 
equation together with the enhanced diagrams of Fig. Hl-b. 



r 



AA/WV 



i 



AAAAA 




/WW 



AAAAA 



Figure 3: The picture of interac- 
tion of two photons with virtual- 
ities larger than a "soft" scale, k 
and k' are transverse momenta of 
quarks. 



Our strategy will be the following: we will solve the non-linear equation assuming that Q' 
and Q ~ l/b t <C k and k' and, using this solution, we will come back to Eq. (j2.1|) and we will 
discuss the contribution of the region where Q + Q' ~ l/b t <C k or /and k'., while Q' ~ Q ~ k 
or k' . 



Finally, Eq. (j2.1J) can be written in the form (see Ref. JH] for details 



dN(k,y;b t 
dy 



(x(m)N(k,y;b t ) -N 2 (k,y;b t 



(2.7) 



where 'y(k) is an operator corresponding to the anomalous dimension of the gluon structure 
function and it is equal to 

m = i + ^ (2.8) 



The form of the first term on the r.h.s. of Eq. ()2.7|) as well as definition of the variable £, 
we will discuss in the next section. Here, we would like only to draw your attention to the fact 
that the emission of the gluons in the BFKL equation is described by the first tern in the r.h.s. 
of Eq. (|2.1|) (see Fig. |2J) which enters at the same value of the impact parameter as the l.h.s. of 
the equation. 

Function ^(7) is an eigen value of the BFKL equation [27] 



X (l) = 2V(1) - ^(1-7) - ^(7) 



(2.9) 



Eq. (|2.7j) we will solve in the semi-classical approach. 



2 These two cases when we can trust the non- linear equation (see Eq. (|2.1(| s | were discussed in the review talk 
of M. Ryskin at DIS'03 conference. We refer to this talk for more details. 
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2.2 The solution to the BFKL equation at large b t 

This solution has been discussed (see Refs. [2HJ HHH EH 1321 EH] ) and at large b t it has a form: 

N(k,k',y;bt) = I 1 p^Mk\r,lH)e &sX< -' r)v - (1 -* /H (2.10) 

J Z7TI 

where £ = hi(bf k 2 k' 2 ) and </>j n is determined by the initial conditions. 

It is easy to see that Eq. ()2.10j) is a solution to the Eq. ()2.1j) without the non-linear term. 
Since the form of solution given by Eq. ()2.10|) is very important both for a derivation of the 
linear part of Eq. (J2.1j) and for understanding the value of the typical dipole sizes in the non- 
linear part of the equation we will discuss this solution in some details follow Refs. [301 EH E2] • 
The general solution to the BFKL equation in the coordinate representation was derived in 
Ref . [321 ) namely (see Fig. EJ), 

N(r u r 2 ;y,b t ) = (2.11) 
/ P~ Mr, ri) d 2 R x d 2 R 2 5(Rx -R 2 - h) e^ y V(r h Rx; 7) V(r 2 , R 2 ; 1 - 7 ) 

with 

w(t) = OLs\{i) (2-12) 

and 



The integration over R\ in Eq. (j2.11jl was performed in Refs. [301 EI] with the result: 
J d 2 R 1 V{n,R x - 1 ) V(rx,\R x -bt\;\-i) = 

( 7 (1 - 7 )) 2 

+ p 2 a: 1 -'V' 1 - 7 F(l - 7, 1 - 7, 2 - 2 7 , x) F(l - 7, 1 - 7, 2 - 2 7 , x*) ) , (2.14) 
where F is hypergeometric function j2H] and x x* is equal to [301 131] 

^,2 ™2 

n* = — U (2.15) 

(6 - z x f ht - z 2 r 2it ) 2 (b - z x r x ,t ~ z 2 r 2 , t ) 2 

and Ci and c 2 are function of 7. 

At large values of b t ^> rx and r 2 the argument x x* is small, namely 

~,2 2 

xx * = 1111 (2.16) 
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It is interesting to notice that xx* turns out to be small and simple also for b t < r 2 with 
7*1 < ?*2. Indeed, for such values of b t 



XX 



-y 2 ,y2 2 
*2 ^2 ' 2 



;2.i7) 



Expanding Eq. (|2.14|) and rewriting the answer in the transverse momentum representation 
given by Eq. ()2.2|) and Eq. (|2.3jl we reduce the general solution of Eq. (|2.14j) to Eq. ()2.10j) . 
Using Eq. (|2.10|) one can see that the linear part of Eq. (|2.7|) indeed gives the BFKL equation. 
Since all derivative in Eq. (J2.7j) are taken with respect to £ we can consider <ft in in Eq. (J2.7j) as 
an arbitrary function of b t . 

The exact form of the initial condition depends on the particular reaction. We choose as an 
instructive example the virtual photon-photon scattering (see Fig- EJ) - 




Fig. Ha Fig. Hb Fig. He 



Figure 4: "Fan" diagrams (Fig. E]-a ), the first enhanced diagram (Fig. 0]-b ) and the first 
correction to "fan" diagrams (Fig. 0] -c ) for 7* — 7* interaction. The thick lines in Fig. 0] -c 
denote the BFKL Pomeron . 

The advantage of this process is the fact that we know the initial condition which is under 
control of pQCD except large bt behavior. The problem of large b t behavior of the Born 
amplitude is addressed in Ref. [22] and it is shown that 0j n can be chosen in the form 

TV" 2 — 1 1 
(j>in(k', y; b t ) = na 2 s ^— (m n b t ) 4 K 4 {2 m w b t ) - (2.18) 

to satisfy both the 1/bf of Born approximation in perturbative QCD, namely, 

/V 2 — 1 1 

<t> in (k',y;b t ) = *c? g -^ - (2.19) 

and the non-perturbative e~ 2m,rbt behavior at large b t ^> l/(2m !r & i ). 



S 



N(r,,r 2 ;y,b) = 

V(r 2 ,R 2 ;l-y) 




CO(Y)V Figure 5: The BFKL 

Pomeron contribution at 
fixed hf. 



2.3 Semi-classical approach 

In semi-classical approach we are looking for a solution in the form: 

N(k, k', y; b t ) = e^ yM) y ~ {1 ~ ^vM) ) € + P{k',h) ( 2 .20) 

where y, b t ) and 7(£, y, b t ) are smooth functions of y and £: du(£, y, b t )/dy <C ^-> 2 (£, y, &t), 
du(£,y,b t )/d£ < w(e,y,6t)(l-7(^yA)), d^,y,b t )/d^ < ( 1 - 7 ) 2 (£, y, &*) and 

dy(£,y,b t )/dy < w(£,y, &t) (1 - 7(6 2/, &t)) 3 - 

Assuming Eq. (|2.2(J|) we can use the method of characteristic (see, for example, Ref. [33] ) 
to solve the non-linear equation. For equation in the form 

F(y,£,S, 1 ,u) = (2.21) 

where S = ujy+j^ + ff, we can introduce the set of characteristic lines on which £(i),y(i), S(t), 
u(t), 7(t) are functions of variable t (artificial time), which satisfy the following equations: 



(4.) ^ = -(F e + 7 F s ), (5.) ^ = -(F y + uF s ), (2.22) 

where F y = dF( -y&Sav) e ^ c 

We can reduce the master equation (see Eq. (|2.1j) ) in semi-classical approach to the form 
of Eq. (j2~2Tj) . namely, 

w - «5X(7) + a 5 e s = 0. (2.23) 

Using Eq. (|2.20j) and Eq. (|2.23|) we can write the set of equations (see Eq. (I2.22jl ) in the 
form: 

-77 = - «s —j , (2.) — = 1 , (3.) — = a s (1 - 7) —\ + w , 

ai a 7 at at 0(7 

(4.) ^ = a s (l- 7 )e 5 , (5.) ^ = - a 5 a;e 5 , (2.24) 
at at 



3 The form of Eq. (|2.2U|) is obtained from the solution of the linear BFKL equation, that has been discussed 
above. 
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We will solve these equations in the next section, but before doing this we are going to illus- 
trate the method of characteristic solving the linear equation neglecting the non-linear tern in 
Eq. (|2.1|) . In semi-classical approach the linear equation has a form: 

u = asxil) • (2.25) 
For linear equation the set of Eq. (|2.24|) reduces to two equations: 

Ty = - as —' ITy = as(1 ~ 7) ^ + W ' (2 ' 26) 
with both 7 and u being constant as function of y. 
The solution is very simple, namely, 

S = x(7s) (V - Vo) - (1 - 75) (£ - £0) + P(b t ) (2-27) 

with 7s given by 

e - £0 = -as ^f^l™ (V ~ Vo) (2-28) 

Comparing Eq. (I2~2TI) and Eq. (I2~2~K1) with the solution to the BFKL equation of Eq. (I2~TU1) 

one can see that the semi-classical approach is the exact solution of the linear (BFKL) equation 
in which the contour integral over 7 is taken by the steepest decent method. Eq. ()2.28|) is the 
equation for the saddle point value of 7 = 75. The accuracy of the semi-classical approach 
is even worse than one for the steepest decent method since we cannot guarantee the pre- 
exponential factor in the semi-classical calculation which appears in the steepest decent method. 
However, the semi-classical approach has a great advantage since it allows us to treat the non- 
linear equation within the same framework as the linear one without major complications. 

To solve Eq. ()2.24|) we need to find out the initial conditions for this set of equations, which 
we derive from the Glauber-Mueller formula for 7* — 7* scattering ( see Ref. [22] ), namely: 

N(x,y = y ,-b t ) = (l- e -^.«W) (2.29) 
where N BA is the dipole amplitude in the Born approximation. It has been calculated in Ref. 

N 2 — 1 1 
N BA (x, y = y 0} b t ) = 7TQ ! |-|^-(m 7r 6 t ) 4 K4(2m 7r 6 t )e« = - f(y = y - b t ) e* , (2.30) 

where the variable £ is defined as 

wffw) • (2 - 31) 

Comparing Eq. (j2.31|) with Eq. (j2.16|) and Eq. (j2.17|) we notice that £ = ln(xx*) with z 2 = 
— z 2 = h. We put these values of z 2 just for the sake of simplicity. 
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One can see that £ in the coordinate representation can be replaced by — £ + 4 In 2 with £ 

£ = In — ^ ^ . (2.32) 



A;' 2 

At large value of bt we have the same definition of £ that has been used but Eq. ()2.32j) gives us 
a generalization which allows us to treat the case of b t = 0. 

Using Eq. ()2.2j) and Eq. ()2.3j) we can find the initial condition in the momentum represen- 
tation: 

iV (£, y = y ; h) = l - T (r(£; bt) = 4 e"« f(y = y ; hj) (2.33) 

where T is incomplete Euler gamma function j2HI of zeroth order. The argument r(£;6 t ) in 
Eq. ()2.33j) can be rewritten as 

Q 2 s (y = y ,k' 2 ;b t ) 



k 2 



T 0,, 



with 



Qi(y = yo, k' 2 - b t ) = 4f r ^ ^ ( k' 2 b 2 t + 1) 



(2.34) 
(2.35) 



bt > 1/k' 



n a 



iN 2 -l 



(2m^b t ) 4 K 4 (2m n bt 



k' 2 blr 0tCr " 3 3N 2 

where we define value of r at y = y on the critical trajectory for £ = £ cr (7 cr ), see 

next section. The behavior of N and 70 — 1 = S% at y = y which follow from Eq. (J2.33|) is 
shown by Fig. 13 . 
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Initial conditions 
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I N o Vo, 
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1.25 ■ 






1 ■ 






0.75 ■ 






0.5 ■ 






0.25 






°0 


0.2 0.4 0.6 0.8 1 



Figure 6: Initial conditions for 
dipole-dipole amplitude No and 
anomalous dimension 70 at fixed 
y = yo as function of r from 
Glauber- Mueller formula for 7* — 
7* scattering j2H] . 



As one will see below we need also to know a ratio 

o x(7o) ~ § 

-DO - : 

1 - 7o 

for finding the solution. This ratio is shown in Fig. [7| 



(2.36) 
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3 Solution at £ > £ sat 

3.1 Fixed QCD coupling 
3.1.1 Solution 



As has been discussed, see P3 El in this kinematic region we can safely use the semi- 
classical approach for sufficiently large values of £. In this section we will find this semi-classical 
solution and will show, that the main qualitative feature of it, is the existence of special, critical 
trajectory which determines the saturation scale. 

From equations 4 and 5 of Eq. (|2.24|) one can see that 

(ko uj , „. 

-r = ~ 3.37 

C17 1 — 7 

Eq. (|3.37J1 has a very transparent physics since it means that the phase and group velocity of 
wave package, defined by Eq. (|2.20J) . are equal. 

Eq. ()3.37|) has an obvious solution 

uj = -a s B (l - 7) (3.38) 

where B is a constant which has to be determined from the initial condition (see Fig. [7| and 
Eq. dHHSJ))- Using Eq. (EOHj) . Eq. (|2~23J) and Eq. (l2~24l - 4 we obtain that 

p- = a s (1 - 7) (S (l " 7) + X(7)) (3-39) 
ay 
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However, before solving Eq. ([3.39)1 . we would like to draw your attention to the fact that 
Eq. ()3.37|) has itself a very interesting solution if we assume that non-linear corrections in 
Eq. (|2.23jl are small but they are valuable in Eq. ([2.24)1 -4 and Eq. (|2.24)) -5. In other words, 
the non-linear term in the master equation is small in comparison with the linear one on this 
special trajectory but the non-linear contributions is essential in the equations for S% and S y 
dependence on y. In this case Eq. (13.37)1 reduces to 

<hda) = x(t) ^ 3 40 ^ 

dj 1 — 7 

The solution to this equation is 7 = 7 cr ~ 0.37 4 . 

The form of the trajectory is clear from Eq. ([2.24)1 -1. namely 

&r = -a s ^p^(y-y ) + (3.41) 

«7cr 

The value of the dipole amplitude iV we can find using Eq. ()2. 24)1 -3 which can be reduced to 

dS q . , . 

Ty=- a * e > (3 ' 42 > 

which has the solution 

N = - °— 3.43) 

N a s (y-y ) + 1 

Since this solution falls down at large y all our assumptions are self-consistent at least at 
as (y — Vo) 3> 1. The trajectory at small values of as (y — yo) should be found from Eq. ()3.39j) . 



3.1.2 Numerical solution 

To find the solution to the master equation we need to solve Eq. ([3.39)1 and obtain j(y) as a 
function of y, substitute this function into Eq. ([2.23)1 and find out N(y) .These two functions 
7(2/) and N(y) will determine the solution on the trajectory that is given by Eq. ([2.24)1 -1. 

One can see from Fig. [Hj that we can divide all possible trajectories in two part with initial 
7o larger or smaller than 7 = 7o, cr - From 7 = 7 cr starts the critical trajectory on which 
l{y) ~^ lar a t large y. The trajectories to the right of the critical one are close to the straight 
lines which are the trajectories of the linear equation. The trajectories to the left of the critical 
line differs significally from the trajectories of the linear equation. In this domain the non- 
linear corrections, induces by partons interaction in the parton cascade, are large and change 
considerably the physics of the QCD evolution traditionally related to the linear evolution 
equation. 

The same physical picture we can see in Fig. El which shows the quite different behaviour of 
the anomalous dimension 7 along trajectories. For the trajectories to the right of the critical 

4 Eq. H3.40JI has been derived in Ref.pQ and has been discussed in detail in Refs. O EE1 EE| • 1 — 7cr is called 
k in Ref.[T], % rU in Ref. [Sj, A in Refs. [S3 EE] and 70 in Ref. |2S] ■ 
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line 7(2/) very rapidly reaches a constant value as it should be for the linear evolution. On the 
critical line 7(2/) approaches 7 cr (see Eq. ()3.4U|) ) but very slowly. For the trajectories to the 
left of the critical line we have a different behaviour and 7(1/) becomes larger than unity which 
is, of course, an indication that we cannot use a semi-classical approach, at least in the present 
form. 

Fig. I1UI shows the behavior of function S(y) (N = e s ) versus y = ln(l/x). One can see that 
the dipole amplitude iV is decreasing slowly on the critical trajectory while it is rapidly falling 
down on the trajectories to the right of the critical line. It is worthwhile mentioning that at 
small y (y = 4- 15) the y- dependence stems mostly from the dependence of j(y) on y ( see 
Fig- El)- We see the manifestation of the slow decrease of Eq. (|3.43jl only at y > 15. 



3.1.3 A general solution 

Since trajectories cannot cross each other we see (i) that the domain to the left of the critical 
line does not know anything about the domain to the right of the critical line; and (ii) that 
the solution to the right of the critical line is very close to the solution of the linear evolution 
equation. Therefore, the semi-classical approach to the non-linear evolution equation leads to 
an ideapQ to solve the problem considering only the linear evolution equation but with initial 
condition not at fixed value of £ = £0 but on the critical line with N is defined by Eq. (|3.43jl . 
For fixed QCD coupling it is easy to find such a solution. 
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Figure 8: Trajectories for 
the non-linear evolution 
equation ( see Eq. ([2.1)1 ) in 
semi-classical approach at 
b t = 0. . 
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Indeed, the general solution to the linear equation has a form of Eq. (|2.1(J|) where 0j n should 
be defined from the initial condition. Using Eq. (|3.43jl and Eq. (13.41)1 we obtain that 

Eq. (|2.1()j) with this <pi n gives the solution to the linear equation with the initial condition on 
the critical line. 



3.1.4 Impact parameter dependence 

Eq. (J2.24)) does not explicitly depend on b t . The entire b t dependence is hidden in initial 
condition (see Eq. (|2.29)l ) both the initial dipole amplitude N Q and the anomalous dimension 
70 depend actually on the variable r as one can see from Eq. ()2.29)1 ). namely, N = iV(r ) and 
70 = T dN(r )/dT . For example, the critical line corresponds to r 0iCr = 0.254 and 70 = 0.248. 

However, when we solve Eq. 1)2.24)1 -1 for the trajectories, we need to know the initial value 
of £ not only of r. Indeed, the equation for the trajectories can be solved as 

£ = -a 3 / dy h fo(ot) • (3-45) 

Jy a J 

It is obvious from Eq. ()2.33)1 that 

Zofa) = \nf(y = y ,b t ) - ln(r 0iCr (y = y )) . (3.46) 

It should be stressed that the b t dependence also enters the definition of £ = In (k 2 k' 2 bf). 
Substituting Eq. ()3.46)1 into the equation for the critical line (see Eq. 1)3.41)1 ) one obtains the 
b t dependence of the saturation scale: 

Q 2 s (y, k' 2 ; b t ) = Q 2 s (y = yo, k' 2 ; b t ) e as <™> (3.47) 

This equation we can rewrite at large b t using the exact expression for Q 2 (y = yo, k' 2 ; b t ) (see 
Eq. (gip ). It has the form 

Q 2 s(v,k ,2 ;b t ) = — -J vra|^-- -(2m n b t ) 4 K i (2mM e as ^^ (3.48) 

k' 2 bf T 0)Cr 3 iV c 2 

In pQCD region Eq. (EPHJl leads to 

^ wkz^^w 1 ^— < "'°' ■ < 3 - 49 > 

As we will see below this expression for Q s leads to the typical radius of interaction which 
increases as a power of energy, namely, 

R 2 (y) oc e l i^ iy - yo) (3.50) 
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in a full agreement with the analysis given in Ref. j2E]. It gives the power-like increase of the 
total dipole cross section, namely 

<T t ot(dipole) = 2nR 2 oc e 5 ~ ( " o) 

in an explicit violation of the Froissart theorem |45j . 

However, the non-perturbative contributions in the initial condition modify Eq. ()3.49j) at 
large values of b t and Eq. ()3.48|) leads to Q 2 oc e~ 2rri7vbt for large b t > 1/2 m^. It is worthwhile 
mentioning that the b t of the saturation scale is the only way how the b t dependence could 
affect the solution in the saturation domain ( to the left of the critical line in Fig. |H] ). We 
will see below that Q 2 S oc e~ 2rri7rbt generates a tot (dipole) oc ln 2 (l/x) in accordance with the 
Froissart theorem |4*o] . 



3.2 Running as 

3.2.1 A general discussion 

In this section we are looking for the semi-classical solution to our master equation (see 
Eq. (|2.1j) ) considering running QCD coupling, namely, as in Eq. (|2.1|) is equal to as = 
&wI7(3? A 2 )) " ^ ne maui uncertainty in taking into account the running QCD coupling is the 
fact that we do not know what is correct way to incorporate running as into the linear BFKL 
equation |3Z1 EEl E23 • However, it was proven, see [T], that we can safely use Eq. (|2.1|) with 
running as in semi-classical approach. 

It turns out that it is useful to search the semi-classical solution in the form 

N(k,k',y;b t ) = a s (Oe S = a s (0 e^' y ' bt)y ~ V-rt'VM)* + , (3.51) 

where u and 7 are smooth function of y and £. 
In this case the master equation has a form: 

w - «5(0x(7) + 4(0 e s = . (3.52) 
Using Eq. (|2.21|) . Eq. (|2.22|) and the explicit form of as, namely 

4 7T 4 7T 

" s(0 = bu - Hk^WW) = &(* - Wx" (3 ' 53) 

we can derive the following set of equation. 



(c , d X (l) . , . dy 

dt ~ R > ~ ~ as{ °—> (2 - } Tt ~ 

at a 7 
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(4.) ^ = -(F 5 + 7 F S ) =«|(0(1 "7)e S - 4(0^x(t) + 4(0 ^ e s ; 

(5.) ^ = -{F y + ^F 5 ) = - 4(0 u;e s ; (3.54) 

3.2.2 Critical trajectory 

There exists a critical trajectory in this case as well as in the case of fixed as- To see it let us 
assume that e s ~ 1 and a$ <C 1. For such 5* we can neglect the non-linear term in Eq. (J3.52j) 
and Eq. (J3.54|) -3 has a form: 

fy =* S (*M + (1-^). (3.55) 

One sees from Eq. ()3.55j) that there is a trajectory with the same equation for the anomalous 
dimension (7 = 7 cr ) as in the case of constant as (see Eq. ()3.4U|) ) on which S is constant. 
Eq. ()3.54j) -4 allows us to determine this constant. Indeed, we can see from Eq. ()3.54j) -4 that 7 
on this trajectory would be constant if [l] 

e S cr = _^X(Hcr)_ (356) 

4tt 1 - 7cr 

The form of this critical trajectory is given by Eq. (j3.54j) -l. namely 

— -j = -«s(0 — ; , ( 3 - 57 ) 

ay d'jcr 

which leads to 



U - i{k'A)f - (Co(6t) - l{k'M)f = I I (y - y ), (3.58) 

b djcr 

where £o(^t) an d yo are determined by the initial conditions. 

Therefore, we have constant S cr and 7 cr on the critical trajectory of Eq. (J3.58|) . One can 
see from Eq. ()3.54|) -5 that u is also approximately constant on this line , namely, 

uj oc e d T«- . (3.59) 
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3.2.3 b t -dependence of the saturation scale 



The saturation scale, which is the value of the typical virtuality k 2 on the critical line, can be 
determined resolving Eq. (|3.58|) with respect to k 2 . It is easy to see that 



QliyM) = A 2 eV(fc(M-f(*',fe)) 2 + s -?\x' cr M\(y-y ) ; (3>60) 

where 

c <h\ ca' h\ 1 (k' 2 mlK 4 (2m w b t )\ 

fo(&t) - €(k,b t ) = In — ► |6t>V2m» - 2m T 6 t , 

for the initial conditions given by Eq. ()2.29j) . However, we will discuss below that we cannot 
solve the non-linear evolution equation with running QCD coupling with this initial condition. 
We need to evolve first the dipole amplitude with the linear BFKL equation to some value 
of y = Din = ln(l/xj n ) and only at y = yi n we can use the expression of Eq. (|2.29|) for the 
initial condition. It leads to substitution N BFKL (x, y = yi n ; b t ) instead of N BA (x, y = yo; b t ) in 
Eq. (|2.29p . For the solution to the BFKL equation we use Eq. (j2.10|) which give us the following 
equation for N BFKL (x, y = y in ; b t ) : 

N BFKL (x,y = y in ;b t ) = ife^^-vo) e ¥ . (3.6I) 

reflects the fact that the BFKL anomalous dimension is equal to \ and lul = ®-sx(\) |2Tj . 
Eq. (j3.61J) leads to the following expression for the saturation scale at y = yi n : 

Q 2 s {y = y in ,k' 2 ] b t ) = (4f) 2 e 2 ^("U /2 (r 0>cr (A ; /2 ^ + l))- 2 (3.62) 

— Ifc>i/k' , /2 ,4 2 ((2m n b t fK 4 (2m n b t )e^^-y^ 2 
k' 2 bfT^ ' 



Icr 



Eq. (|3.62j) translates in 



r -\n(bfk' 2 A 2 ) for (1/2 m,) > 6 t > l/k 1 
HQ 2 s (y = y m ,k' 2 ;b t )/A 2 ) = (3.63) 

[ -4m 7r 6 t for 6 t > (1/2 m^) 

Eq. ()3.60|) shows that the saturation scale does not depend on b t for b t < b (y) with b 
which is the solution of the following equation 



HQ 2 s (y = y in ,k' 2 ;b )/k 2 ) = ^\ X ' 7 J (7cr) (V ~ yo). (3.64) 

From Eq. ()3.63|) one can see that bn oc exp (hJ^- \x' lcr (lcr)\ (y ~ yo)) at 1/2 ^> bo ^> l/k' 

and b = (l/4m*) V^jx^7c^)l (y - yo) for 6 > 1/277V In other words, Q 2 -> 1/6| for 
6 t > bo if l/k' < b < 1/2 and the saturation scale Q 2 —> e ~ im ^ b t f or 5 Q > \/2m^. 
Fig. ^2 illustrates such a behaviour of the saturation scale. It should be stressed that we have 
a different large b t behaviour for the running a$ due to different initial conditions. 
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Figure 11: The behavior of 
the saturation scale as func- 
tion of y at different val- 
ues of the impact parame- 
ters ib t ). As one can see at 
large values of y the satura- 
tion scale does not depend 
on bt in a striking difference 
with the case of frozen QCD 
coupling. 



3.2.4 Numerical solution 

Our above discussion we support by numerical solution of Eq. ()3.54|) . We simplify a bit this 
set of equations neglecting a| contribution in Eq. (J3.52JI and a| term in Eq. (|3.54|) -4. Fig. IT21 
shows that the trajectories for the running QCD coupling have the same pattern as for frozen 
a s : 

• They can be divided in two parts by the critical line which has the form of Eq. ()3.58|) 
starting from £ > 5. The anomalous dimension on this line rapidly reaches the value of 
7 cr ( see Eq. (|3.4U|) and Fig. [13] ) while S(y) is constant at large y (see Fig. [Hj): 

• All trajectories to the left of critical line in Fig. ED go to the negative values of £ with 
unreasonable values of the anomalous dimension j(y) — > 1.62 (see Fig. IT3*j) and with very 
fast increase of e s ^ shown in Fig. El Such a behaviour reflects the fact that we cannot 
apply the semi-classical method for the solution inside of the saturation (CGC) region 
at least when the trajectories go far away from the critical line. Comparing Fig. |§] and 
Fig. |S1 with Fig. EI] an d Fig. ^] we conclude that solutions for both frozen and running 
QCD coupling have the same qualitative features; 

• Among all trajectories to the right of the critical line we can conditionally separate 
two groups. The trajectories of the first group envelope the critical line and only at 
some value of y they go apart. For the large values of y > y cr they coincide with the 
trajectories of the linear evolution equation. Therefore we can find the dipole amplitude 
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on such trajectories just solving the linear evolution equation with the initial condition 
N = a s e Scr at y = y cr with e Scr given by Eq. ()3.56|) . The trajectories which do not touch 
the critical line belong to the second group. They are the same as trajectories for the 
linear equation. 

• The trajectory which separate these two groups of trajectories to the right of critical 
line is denoted by dotted line in Fig. El and it touches the critical line in the only one 
point. The equation of this tangent line was found in Ref.pQ and it reads: 

y - m '^U^itrJ^' (3 ' 65) 

where y(£, y) and y)) are the coordinates of the point where the tangent line ( dotted 
line in Fig. 112)1 . passing through the point (£,?/), touches the critical line. 



Trajectories 




Figure 12: Trajectories 
for the non-linear evolution 
equation ( see Eq. (|3.54)l ) in 
semi-classical approach at 
b t = for running QCD 
coupling. he dotted line cor- 
responds to the trajectory 
which is tangent to the crit- 
ical line. 



The essential difference between frozen and running QCD coupling occurs in the initial 
condition for the non-linear equation. For the frozen as we used the initial conditions given 
by Eq. 1)2.29)1 and Eq. ()2.3U|) . However, we found out the the non- linear equation with running 
as cannot be solved with the same initial condition. It turns out that we need to evolve our 
parton cascade using the BFKL equation (see Eq. ()2.10p ) for y less than y (x > Xq) and 
use Eq. (12.29)1 as the initial condition at y — yo (x — xq) with N BFKL instead N BA . In our 
numerical calculation we use Eq. (J2.10)) and Eq. ()2.18j) with xq = 1CT 1 to obtain the solution 
that has been discussed. 
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Figure 13: The anoma- 
lous dimension 7(7/) on the 
trajectories for the non- 
linear evolution equation 
(see Eq. The dotted 

line corresponds to 7 = 7 cr . 




y 
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4 Semi-classical Approach in Saturation (Colour Glass 
Condensate ) Domain 



4.1 Reduction of the non-linear equation to the semi-classical form 

As have been discussed we cannot expect to solve the master equation in the saturation (CGC) 
domain considering the dipole amplitude N = e s such as S being smooth functions of y and 
£. Therefore, the fist question, which we have to answer, is for what function we are going to 
use the semi-classical approximation. We follow the idea of Ref. ^H] and introduce the new 
function $ as 

N(y,£;bt) = i [ d? U - e"*^*)) . (4.66) 



We consider $ as a smooth function such that $^ <C $ yy «C $y$y as wen as 

$^ < <&e$j,. Therefore, 

' e -* = (-$ f ) n e"* . (4.67) 



(90 

Substituting Eq. ()4.67|) in Eq. ()2.7|) we reduce it to the following equation for function $ 



^^, (l -^,- ) + 5 ,^^, (4 , 8) 



where L(f) = — /) — 1)//- In double log approximation of pQCD which was discussed 

in Ref.fTT?] function L is equal to zero since %(1 ~~ f) — !//• F° r the derivation of Eq. ()4.68|) 
one needs to take into account that the series for L(f) starts with f 2 term. 

In this section we will solve Eq. (I4.68|) using the powerful idea of the geometrical scaling 
behaviour which will allow us to simplify this equation. 



4.2 Geometrical scaling 

Inside of the saturation domain we expect so called geometrical scaling behaviour of the dipole 
amplitude, namely, this amplitude depends on one variable r = Q 2 s (x; b t )/k 2 instead of three: 
x, k 2 and b t . The physics of such behaviour has been discussed for long time starting from the 
first papers where the saturation scale was introduced [U 121 E] • Recently, it has been found that 
the HERA data supports this idea in the wide range of low x ( x < 1CT 1 ) (see Ref. (HIl for 
details) as well as the numerical solutions of the non-linear evolution equation show this scaling 
behaviour [221 12H HI] ■ However, the theoretical proof I2D1 HHj of this scaling behaviour is 
still in a very embryonic state. The only result, that has been proven so far, is that in the 
semi-classical approach the master non-linear equation has a geometrical scale solution both in 
the saturation (CGC) region jM] as well as for [201 P 5 (Qt/A 2 ) > Q 2 > Ql 

5 It is interesting to mention that the semi-classical solution in Ref. pp shows the geometrical scale behaviour, 
but it was found out only after an explanation given in Ref. |2(J| . 
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In this paper we will find the semi-classical solution to Eq. (J2.7|) inside of the saturation 
(CGC) region. 

We start with the general arguments why the semi-classical leads to the geometrical scale 
behaviour of the dipole amplitude given in Ref. [3*4*] . 

The solution to the linear equation can be written through the moments of the dipole 
amplitude, namely, 

N(k,k',y;b t ) = [ ^0(tf |W ;&^»-< 1 -'rW)« (4.69) 

where 7(0;) is a solution to the following equation: 

a s x(7H) = oj. (4.70) 

In to representation Eq. ()2.7|) has a form 

(u - a s x(l(u))) N(u,&bt), = - a s [ N \J , £; b t ) N (u - c/, £;&,). (4.71) 

J 2 m 

In the semi-classical approach we can use the steepest decent method to take the integral 
over us' in Eq. (|4.71|) and the saddle point in this integration is equal to u>' s = uj/2. After doing 
this we obtain: 

(u - a5X(TH))#>;^)^ (1 - 7MK = (4-72) 



<xs<f>(k',—',b. 



w . , 1 



-2(l-7(f))? 



2' lJ 2V27r 7i ^f 
In Eq. (|4.72|) one can see two different regions in to: 

1. uj ^> Ucrit- Here we can neglect the r.h.s. in Eq. (|4.72jl and anomalous dimension 
j(uj) can be calculated from Eq. ()4.70j) . 

2. uj <C Merit- The r.h.s. of Eq. (|4.72jl is not small and at the fist requirement for the 
solution we need to have 

(1 - >y(u)\ u<Ucra ) = 2(1 - 7(|)U«^J • (4.73) 
This equation has a simple solution jHlj 

(l-T(w)U w J = Cu. (4.74) 

3. Both solutions should match at a certain value u = uj cr u. Therefore, we have 

(1 - 7H)| w>Worii = (1 - 7M)U< Wcr , ( ; (4.75) 

duj ^ UJcrit dijj ^ UJcrit ' \ IG ) 

ttsX(7cr) 

One can see that Eq. (|4.77J) is the same as Eq. ()3.40j) with j(uj cr it) = 7 cr . 
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Of course, we can see that (1 — 7) = Cuj directly from Eq. (|2.24j) but we consider the above dis- 
cussion (see Ref. [33] for more details) as more general since we use the very general properties 
of the semi-classical approach but not the method of characteristics. 

The fact that (1 — 7 oc u in the saturation(CGC) domain means that the dipole amplitude 
is actually a function of one variable 

z = a s * (7cr) s V ~ t + P(k',b t ) . (4.79) 

(1 - lev) 

One can derive this directly from Eq. (J4.69j) substituting Eq. (I4.78j) and changing the inte- 
gration variable from ui to uj' = ocsxilcr) UJ - 

4.3 Solution 

Looking for the function of the new variable z given by Eq. (|4.79|) we reduce Eq. ()4.66|) and 
Eq. (jOHl) to the form: 

N(z) = ~ dz' (l - e-* (2,) ) ; (4.80) 
2 Jz 



We can rewrite Eq. (j4.81|) as a set of equations: 



dz 



D(z) ; (4.82) 



VI - 7cr cLD($) / ^ ' d® ' 1 ^ 

which can be easily solved. It should be stressed that Eq. ()4.83j) does not depend on parameters 
and, therefore, it could be solved numerically without loosing any generality. 

Integrating Eq. (J4.83|) over D and $ and introducing a new variable 7 cr = 1 — 7 cr we obtain: 

v/™ \ D 2 — D 2 rD n 

^i^—J^L + x(D) D - X (lcr)%r - I x{D')dD> + In if = (4.84) 

7cr ^ J %r M) 

$ + e- $ - $ - e-*° , 
where -Do = D(zq) and $0 — ^(-^o) are initial values for -D(z) and $(2). 

We can see two distinct regions of <£> where we can obtain an analytic solutions to Eq. (|4.82j) 
and Eq. fEEfl . 
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z»z where z is the value of z for the critical line as has been discussed. In this 
case D is approaching the singularity D = 1 in the r.h.s. of Eq. (|4.84|) . Therefore this 
equation is reduced to 

' : $ , (4.85) 



1 - D 

or D = 1 - 1/$. Solving Eq. (J4~%2|) we obtain $ = z - \nz. Therefore, Eq. (|P5j) 
leads to the following dipole amplitude for large z 

N{z) = \ f (l - z'e- z ') dz' . (4.86) 

Therefore, we see that at large values of z our function <l> is also large and we can safely 
use the semi-classical approach. 

• z — > zo . For such z we cannot use the semi-classical approach for function $. As we 
have discussed we can use here rather the semi-classical approach for the dipole amplitude 
N, which is equal to 

N(z)\ z ^ zo = N(z )e^ z (4.87) 

near to the critical line. 

However, if <3> = $(zo) turns out to be so small that we can expand Eq. (|4.66|) we can 
obtain from Eq. (I4~%7I) that 

$(z)U* = 2% r N(z )e %rZ (4.88) 
In this region we can re-derive Eq. ()4.68|) replacing $^ by (ln$)^. 

In Fig. [T^lwe presented the numerical solution to Eq. ()4.82|) and Eq. ()4.83|) . We use N = 
0.32 as the initial condition which comes from the solution to the right of the critical line for 
running QCD coupling. 

One can see that the exact solution can be well approximated by the asymptotic formula 
$ = z — ln(z) starting from z — 1. It should be stressed that this formula leads to a quite 
different asymptotic behavior than it occurs in the so called double log approximation of pQCD 
(DLA) (see Ref. ^H])- In DLA it turns out that $ oc z 2 deeply in the saturation (CGC) region. 
For lower value of z we can use Eq. (J4.88|) as seen from Fig. • 



4.4 Saturation for running a s 

As has been discussed it is not theoretically clear how to include the running QCD coupling in 
the non-linear evolution equation (see Eq. (|2.1)) ). In Ref. pQ it was shown that we can safely 
consider as a function of the initial dipole Xqi for the kinematic region to the right of the 
critical line. However inside the saturation region we certainly cannot do this. We decided to 
consider the running as in Eq. (j2.1j) to be frozen on the critical line 6 . We cannot prove this 

6 In other words we take as = &s{Q 2 s {y)) inside of the saturation region. 
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Figure 15: The semi- 
classical solution to the 
non-linear equation as a 
function of the geometri- 
cal scaling variable z = 
ln(Q 2 s (x;b t )/k 2 ). 



assumption but it looks natural from the point of view that physics of saturation is determined 
by one scale: the saturation momentum. 

In framework of this approach we obtain the same set of equations (see Eq. ()4.82|) and 
Eq. (|4.83jl ) as for fixed as. The only difference is the new definition of the variable z. Instead 
of Eq. (gZZSD we have 



X{lc 



V 1 - lev) 

Using Eq. ()3.60|) we can rewrite Eq. ()4.89|) in the form: 

z = \n(Q 2 s (y;b t )/k 2 



y - e + P(k',b t ) . 



(4.89) 



(4.90) 



However, the statement that Eq. (j4.90J) as well as the fact that we have the same set of equations, 
is correct only for high energy ( large values of y ) for which 



—\x', /cr (icr)\ (y 



vo) » (£o(M - Z(k',b t )f 



5 Unitarity Bound 

The solution of the previous section allows us to calculate a high energy behavior of the ob- 
servables. 
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xG(Q 2 ,x). The dipole amplitude in momentum representation, that has been calculated, 
is closely related to the gluon structure function, namely: 



xG{Q 2 ,x) 



2W « fQ \e ' - 



as it Jo 



d%N{y,Z,b t ) ■ 



(5.91; 



The simplest way to understand Eq. (j5.91j) is to notice that at short distances the dipole 
amplitude in the coordinate space is directly related to the gluon structure function. Indeed, 
the dipole total cross section is equal to 



crdipoie{x,r t ) = — ^-^vr r t xG{ — ,x) 



2N(x,y;b t )=2r 2 



kdk f d 2 b t J {kr t )N{k,y;b t ) . 



(5.92) 



The last relation is Eq. (J2.2|) in which we use a new notation for the dipole size (r t ). Using 
moment representation for both xG and N in Eq. (|5.92j) . namely 



xG(Q 2 ,x) 



2ni 



(5.93) 



we can find the relation of Eq. (J5.91|) repeating simple calculation given in Ref . [13] ■ 
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Figure 16: The ratio of the 
gluon structure function in 
the saturation (CGC) re- 
gion to it's asymptotic value 
given by Eq. (|5.95|) as a 
function of energy (y = 
ln(l/x)) at different values 
ofQ 2 . ' 
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Using asymptotic expression for N = (see Fig. Ho")) one can calculate xG asymp using 
Eq. (j5.91|) . We can take all integrals analytically if we assume a simplified version of Eq. ()3.48|) . 
namely, Q 2 (y; b t ) = Q 2 (y; b t = 0) Exp(-2m^b t ) or 

z = as ^ la) V ~ In k 2 - 2 mA ~ Zq = \n(Q 2 s {y) / k 2 ) - 2mA - z . (5.94) 

Integrating over b t in the saturation region [z > or Q 2 < Q 2 (y, b t = 0)) we obtain 

7ro;5 4m* 

Fig. [THl shows the ratio of xG calculated directly from Eq. ()5.9H) using the numerical solu- 
tion which we have discussed in the previous section to the asymptotic expression of Eq. ()5.95|) . 
Actually only to find out this difference we developed the method of solving for the master equa- 
tion since the asymptotic behavior is clear from the general properties of the dipole scattering 
amplitude, namely, from the fact that N(r 2 ,y,bt) < 1 due to the unitarity constraint. 

In Eq. (j5.95|) we consider so high energy that the upper limit in bt integration which is equal 
to bt = &o — jm~ m (Qs(?/)/^ 2 ) i s larger than l/2m 7r . In the region where bo < l/2m 7r we can 
use the pure pQCD expression for the saturation momentum given by Eq. ()3.49|) . The value of 
&o is equal to 

J.4 1 2 ( N c - 1) a s p^- (y-y ) / r r\p \ 

b ° = k^ na s 6N 2 e s ^-^yy . (5.96) 



c 



Using Eq. ()5.96|) and Eq. (j5.91j) we have 

xG asymp (Q 2 ,x) = (5.97) 



( N c - 1 ) Q c a 3 ± ^(y-yo) 



J- f Q dk 2 2[ b ° db 2 Hb 2 /b 2 ) = 4 
iias Jo Jo \ 6 N£ 7r k 

Therefore, the gluon structure function shows a power-like increase for y smaller y a given 
by equation b (y a ) = 1/2 m^ with b taken from Eq. ()5.96|) . 

For running as we have to resolve Eq. ()3.64|) to find b (y) and substitute it to Eq. ()5.97|) . 
It is easy to see that in the pQCD region 

^2 _ \ f s e Wl ( yin ~ yo ^ e?J nr x '~fcr (tct) {y - yo ) _ (5.98) 

However, at y = y a bo from Eq. ()5.98|) reaches the value bo = an d for y > y a the value of 
bo is equal to 

bo = ^x'.Jicr) (y-yo). (5.99) 
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cr(dipole). This cross section is equal to 



a(dipole) 



d 2 b t N{r t ,y;b t 



(5.100) 



We use Eq. (12. 2 J) to calculate N(r t , y\ b t ) and it is easy to do for the asymptotic formula for 
N(k,y;b t ) — \z = ln(Q s (y;b t )/k) = K (k/Q s (y,b t )). Here, we replace log contribution by 
McDonald function which reproduces this log at small value of the argument and leads to an 
exponential decrease at A; > Q s - For such large k we cannot trust our simple formula and 
have to replace it by the perturbative QCD expression which will give a small contribution in 
the total cross section due to the fast decrease at large bt. The integral over k in Eq. ()2.2|) 
give N(r t ,y;b t ) = 1. Integration over b t we can take using z from Eq. (J5.94j) with changing 
k — > 2/r t . Therefore the asymptotic behavior for a(dipole) is 



aidipolt) 



2tt 
Ami 



\n 2 (r 2 Q 2 s (y,b t = 0)) . 



(5.101) 



For y < y a the upper limit in b t integration should be taken from Eq. ()5.96|) and 



a(dipole) = 2nbl(y) 



-7T r 1 r 2 a s ^ 



7T- 



\ as 



(y-vo) 



6N3 



(5.102) 



7*— 7* cross section. 

The total cross section for 7* — 7* interaction can be written as follows (see Fig. |2J): 

a(Q 1 ,Q 2 ,W) = (5.103) 

J dzi J q dz 2 J d 2 r 2 u J d 2 r 2tt \^(Qi, z 1 ,r ltt )\ 2 \^(Q 2l z 2 ,r 2jt )\ 2 a(dipole) 

where aidipole) is given by Eq. ()5.100|) while the wave functions of virtual photons are well 
known 



The asymptotic expression using Eq. ()5.101|) has been written and discussed in Ref. 
Here, we can add to discussions in this paper the fact that the asymptotic behavior starts rather 
early as one can see from Fig. namely, at y > 2 4- 2.5. 

All above estimates were done for the frozen QCD coupling because only in this case we 
obtained the semi-classical solution in entire kinematic region. For running as we have not 
solved the master equation since (i) we could not prove the geometrical scaling behavior in the 
saturation (CGC) domain; and (ii) we do not know what is a correct way to take into account 
the running QCD coupling for the BFKL emission. However, all asymptotic formulae have so 
transparent physical meaning that we can write them for running QCD coupling as well. 

Indeed, Eq. (I5.95|) has three logs: 

• One comes actually from the fact that the dipole amplitude in the coordinate space 
should approach unity at high energies. This is a direct consequence of the unitarity 
constraint and, because of this, will be the same for any solution to the master equation; 
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• Two logs are originated from the impact parameter integration. The limit for this 
integration stems from the kinematics of the saturation (CGC) domain, namely, from the 
fact that k 2 < Q 2 (y,b t ) in this region. This relation leads to the b t < j^r m (QsA 2 ) as 
we found in Eq. (|3.60j) . 

Therefore, we expect the same form of the asymptotic behaviors for running «5 as for frozen 
QCD coupling in terms of the saturation scale and the size of scattering dipoles. However, the 
energy behavior of the unitarity bounds will be quite different since the saturation scale is 
proportional to Q 2 oc ^/ln(l/x) for the running QCD coupling while it increases as ln(l/x) for 
constant ag. 



6 Accuracy of the approach 

As we have discussed we neglected the contributions to the non-linear term in Eq. (j2.1j) from 
the region of integration xu ~ £02 ~ 26 t 3> To estimate the accuracy of our approach we 
have to calculate the term that we neglected using the solution to the non-linear equation that 
we found in this paper. 

We need to calculate the diagram of Fig. [TJI which is equal (see Eq. ()2.6|) and discussion 
around this equation ) 

2 

a s (x 01 ) J d 2 Ab t * 01 N 2 (2(b t + Ab t ),y,Ab t ) . (6.104) 

For b t ^> bo(y) where b${y) is the radius of interaction that stems from our solution, it is easy 
to see that Eq. (|6.1U4j) leads to the contribution which is of the order of 

since the dipole amplitude N is of the order of unity for b t < b . One can see that AiV <C N 
and we can safely neglect this contribution. 

For r 2 < b t < b (y) Eq. (16.10411 gives 

AiV « a s (x 01 ) -gL . 

Since X01 bo(y) and bo(y) increases at least logarithmically with energy AiV is very small. 
To the right of the critical line the non-linear term of the evolution equation is small and 

2 2 

the diagram of Fig. El has an additional smallness which is proportional to ce six 01) x °}[ 2 for 
b t > r 2 . 

It is interesting to notice that we can repeat the same estimates for the linear term in 
Eq. (|2.1J) replacing N 2 in Eq. (|6.104jl by 2N. These estimates would lead to the same result. It 
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X 01 - r i 



N'(2b„y,Ab t ) 



Figure 17: Corrections to 
the non-linear equation due 
to interaction of two large 
dipoles (xi2 ~ X02 ~ 
2b t > r 2 ). 



means that the large bt dependence of the linear term (see, for example, Eq. ()3.49|) ) stems from 
the dipole sizes smaller than the size of the target. This is a transparent explanation why the 
first linear term in Eq. (j2.1)l has the form given by Eq. ()2.7)1 without any additional assumption 
about the size of the dipoles in the non-linear equation 7 . 

The most important restriction of the applicability of Eq. (J2.1)) stems from enhanced diagram 
of Fig. I3J-C. Estimating these enhanced diagrams we should distinguish two different kinematic 
regions : k 2 > Q 2 and k 2 < Q 2 . In the first one the integration over 1/2 in Fig. HJ-b converges 
at 1/2 oc 1/ds "C y and, therefore, the enhanced diagram reduces to the diagrams of Fig. HJ-a ( 
see Refs. [U El HH1 EEl Ej for details). In this region the nuclear target has an advantage that 
the 'fan' diagram originated from the enhanced diagram has an additional suppression with 
respect to parameter as A3 [T7| 13%]. 

However, the most interesting region is the second one (k 2 < Q 2 )- In this region there is 
no reason to expect that the integral over y 2 will be concentrated at y2 <C y and the enhanced 
diagrams will lead to the dominant contribution. The nuclear target has no advantages in this 
region since the solution here does not depend on the initial conditions. As one can see directly 
from Eq. (|2.1|) the particular properties of target enters only in the initial conditions. We are 
concentrated our efforts on dealing with this kinematic region in this section. 

The expression for the first correction due to enhanced diagram to "fan" diagrams presented 
in Fig. El-c is the following [TJ [T7] : 

J d%AN enhance = H| J d%£ d yi J d 2 k 1 d 2 b l dy 2 J d 2 k 2 d 2 b 2 (6.105) 

N(k, k u y-yt, b^N 2 ^, k 2 , y x - y 2 , b 2 )N(k 1 , k', y 2 , b 3 ) 

All notations are clear from Fig. EJ-b and Fig. 0J-c. It should be stressed that all amplitudes 
in Eq. (J6.KJ5)) except the last one(A r (fc 1 , k', 1/2, 63)) are different from the solution of the non- 
linear equation (see Fig. EJ-c). They actually satisfy a different equation which it is easy to 

7 We thank our referee who pointed out this transparent illustration of the form of Eq. (|2.7J) . 
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write. Indeed, the slight glance at Fig. HJ-c shows that the equation for N(k, ki, y — yi, &i) can 
be written in the form: 

dN(k,k u y -yi]b t ) 



dy 

&s ( x(j(k))N(k, ki, y - yi] b t ) - N{k, h,y- y x ; b t ) N{k, k', y; b t ] 



(6.106) 



where N(k, k', y\ b t ) is a solution of Eq. ()2.7|) . Since the solution to this equation is a function 
of one variable z and at y± — it should coincide with the solution to Eq. (|2.7j) one can find 
that N(k,k 1: y — yi,b t ) = N(z — Z\) while N(k, k' ,y;b t ) = N(z) where N(z) is defined in 
Eq. (|4.80|) . We can rewrite the argument of N(k, ki, y — y\ \ bt) in more convenient form using 
Eq. (|3~MJ) and Eq. (loHHj) . namely 



'F( y- yi )S (b t ) 
k 2 k\ b\ 



N(k, kl , y- yi -b t ) = N[ '^Ji^ (6-107) 



where F(y — yi) is a function of rapidity. This function is different for frozen and running 
as, namely F{y) oc e XrV for constant as and F(y) oc e Arv ^ for running QCD coupling (see 
Eq. ()3.64|) and Eq. ()5.98|0 . S(b t ) is a non-perturbative profile function for which we assume that 
S(b t ) — > e ~ 2m T b t a ^ large b t {b t > 1/2/71^). Using Eq. ()6.107|) we can evaluate all integrals in 
Eq. ()6.105|) . One can see that in pQCD region for not very high energies where we can consider 
that S(b t ) = 1 we have 

J d%AN enhance cx ^ ln2( ^ /2) J q V d Vl J q V1 dy 2 F{y - Vl ) F( Vl - y 2 ) F{y 2 ) (6.108) 

Calculating the ration of AiV to iV to characterize the scale of corrections we obtain: 

fd%AN g| d Vl S» 1 dy 2 F(y - Vl ) F( yi - y 2 ) F(y 2 ) 
Sd*b t N K N* F(y) ■ (b ' iUJj 

For frozen as this ratio is of the order of jjiy 2 8 which means that we can neglect the 
enhanced diagrams only in limited range of energies. However for running as this ratio is 
proportional to j^e^^^ 1 ^ and the range of energy (rapidity) where we can still neglect the 
enhanced diagrams shrinks even more than for constant as- 

At first sight the above discussion is in the direct contradiction with the paper of Mueller 
and Patel [17] where was shown that (i) the enhanced diagram of Fig. HJ-b has no extra N c 
suppression with respect to 'fan' diagrams of Fig. HJ-a; and (ii) the normalization of iV could be 
chosen in a such way that extra factor 1/N 2 will not appear in the calculation. 

To clarify our calculations we will repeat them in a toy-model: 'fan' diagrams for Pomeron 
exchanges with a' P = and intercept «p(0) = 1 + A. For the nuclear target such model 
correspond to Schwimmer resumation and describes the hadron - nucleus interaction at 

8 We do not consider ln(fc 2 /fc' 2 ) here as a large factor, since in the saturation region this factor should not 
be specially large. 
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high energy. In this model the dipole-dipole interaction is given by sum of 'fan' diagrams of 
Fig. HJa, namely, 



9\ 92e Ay j/>i giA 

1 + ^,^(6^-1) G 3p 



m = - , : ":: 2 ; XM ^ ^ ^ = i (6.110) 



where g\ and g 2 are vertices of Pomeron interaction with dipoles 1 and 2 [r\ <C r 2 ). Actually, 
gi contains an extra factor 1/A" 2 as one can see from Eq. (j2.18j) . In Eq. ()6.110j) we assumed 
that N(y) — > 1 at y — ► oo as it follows from Eq. (|2.1|) . The subset of 'fan' diagrams that starts 
from the Pomeron and finishes at dipole 2 (target) leads to the answer 

*tf-°> = i + ^(^°- 1) ^ Ae(!, '- 0) ' (6 ' U1) 

where O(y) is equal to 1 for y > being zero for negative y. 

One can see that N(yi — 2/2) starts with one Pomeron exchange e A ( yi ~ m \ Summing all 'fan' 
diagrams attached to this Pomeron we have 

N(Vi-V2) = e A( « } £ ( A (2/1-2/2))" = 1 ; (6.H2) 

n=0 n - 

while 

N( y - Vl ) = gi e A ^ (&(y-yi)T = 9i ; (6.H3) 

n=0 n - 

N(V2~0) = r ^ A. ( 6 .114) 

Collecting Eq. (jfi.lllj) - Eq. ()6.114j) we obtain for Eq. (j6.105|) and taking into account that 
77^- = 1 we obtain 



and 



AN' 



enhance _ ^3P f . f , « 2 



I cfoj I dy 2 oc^y 2 (6.115) 



which coincide with Eq. (|6.109|) . 

Formally speaking the enhanced diagrams is suppressed by factor 1/N% and have not been 
included in the non-linear equation as well as a large number of other diagrams of the order of 
1/Af?. Therefore, our estimates show that 1/N C corrections could be rather large. Actually, we 
know how to calculate such kind of corrections ( see paper of Balitsky in Ref. ^7j and paper 
of Wiegert in Ref. [To]) 9 

9 We thank Heribert Weigert for discussion of 1/N C corrections. He shared with us his first attempts to 
calculate numerically these corrections. 



34 



7 Discussion and Conclusions 



In this paper we discussed the semi-classical approach to the solution of the non-linear equation 
which has been suggested and developed before [H HU HH1 HHj- The new domain which we 
explored in the semi-classical approach is the saturation (CGC) region which was out of reach 
by this method before. The new result is also related to the impact parameter dependence 
which has not been considered in Refs. El G2] and has been approach only in the simplified 
case of double log approximation of pQCD in Ref. [T9] . 

We hope that we demonstrated the main advantages of the semi-classical approach, namely 

• The natural appearance of the new saturation scale as the momentum of the critical 
line which divides the set of trajectories in two parts: the first group is trajectories which 
approach the trajectories of the linear equation at large value of £ = \n{k 2 k l2 bf) (see Fig. El 
for notations); and the second group is trajectories that go apart from the trajectories of 
the linear equation; 

• The solution in the saturation (CGC) domain reproduces the geometrical scaling be- 
havior , namely, cr dipo i e (r t , x; b t ) turns out to be a function of one variable r = r 2 Q 2 (x; b t ); 

• The impact parameter dependence of the saturation scale stems from the initial condi- 
tion (Born approximation) as has been discussed in Refs. [221 121] but not from the kernel 
of the linear evolution equation (BFKL kernel)) as was suggested in Ref.(2S]. The large b t 
behavior of the saturation scale should be determined by non-perturbative contribution 
to the Born amplitude which has been considered in Ref. [25]. 

We found asymptotic behavior both for the dipole amplitude N and function $ (see Eq. (|4.66)l ). 
It should be stressed that the behavior of function <3> at the large values of the argument z is 
quite different from the double log approximation case discussed in Ref.jTH]. It should be 
stressed that we found that the exact solution could be approached by asymptotic one starting 
with y = ln(l/x) rs 2 4- 2.5. 

We also found the unitarity bound for the gluon structure function and for the dipole cross 
section, given by Eq. (|5.95|) and Eq. (|5.1U1|) . which depends on the energy dependence of the 
saturation scale. The fact that high energy behavior of the dipole cross section depends on 
In (Q 2 (x; b t = 0) rf) differs our analysis from the considerations discussed before [221 12H I2H] ■ It 
should be stressed that our formulae lead to quite different dependence on energy for the high 
energy asymptotic behavior. Even for the frozen QCD coupling our results are different from 
the estimates in Refs. (221 [21] • For example, we predict that 




(7.116) 



instead of 




(7.117) 
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Figure 18: The dipole ampli- 
tude as a function of the im- 
pact parameter b as it emerges 
from discussion in section 6 (see 
Eq. (IFHU3j) ). The dotted line 
shows the correction term evalu- 
ated using Eq. (jbMU4|) . 



as was obtained in Refs. [23| 124] . It is worthwhile mentioning that 7 cr ~ 0.37 (see Eq. ()3.40|) ) 
is not related to the running QCD coupling. It stems from the full BFKL kernel while 7 cr = | 
arises only if we approximate the full BFKL kernel by the double log formula x{l) = 1/7- 

In the case of running as the difference is even more remarkable since our formulae lead to 
a different energy behavior. Indeed, we have 

a(dipole) < _^ p£ 2^ (7 . 118) 

as one can see from Eq. f)3.60j) . Notice that extra factor 4 in the dominator stems from the 
different initial conditions that we used in the case of running as- 

Actually, the discussion about corrections to the non-linear master equation given in section 
6 can be summarize in the picture shown in Fig. The corrections from the dipole which 
sizes larger than the size of the target lead to a large b tail but it starts at very low value and 
it could not contribute to the energy behaviour of the total cross section. The entire picture is 
the same in the main features as it was suggested in Ref. j21] 

We started this paper having in mind two questions to answer. The practical one: to find 
the analytical solution which everyone can check and which includes the impact parameter 
dependence. We firmly believe that only after finding such a solution we could develop the 
numerical methods of solving our master non-linear equation (see Eq. ()2.1Jl ). Since the semi- 
classical approach is valid for the large values of y = ln(l/a;) ^> InQ 2 we can use this solution 
as a guiding one for the numerical attempts to solve the non-linear equation. 

Our second motivation was to reach a theoretical understand how to incorporate the non- 
perturbative QCD corrections in the nonlinear equation. These corrections have to be essential 
at large values of the impact parameters where they will change the power-like decrease in pQCD 
to the exponential one as follows from the spectrum of the observed hadron j^SJ. We believe 
that we answered this question by demonstrating that we should take the non-perturbative 
correction only in initial conditions (in Born amplitude). How to take into account these non- 
perturbative corrections in the Born amplitude was considered in our previous publications 
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How to incorporate the non-perturbative corrections in initial conditions or in the kernel 
actually crucially depends on the strategy of the approach. If you consider first the solution 
of the linear equation you have to include the non-perturbative corrections at large b t in the 
BFKL kernel as was discussed in Ref. j2Hj. We believe that the correct strategy is to solve the 
non-linear equation. As has been discussed for long time PJ EHJ EE1 IS] the non-linear dynamics 
leads to a suppression at low x of the density of partons with low transverse momenta (pt). 
Since due to the uncertainty principle Ab t pt ~ 1 only partons in the beginning of evolution (at 
x ~ 1) could have large values of b t . In other words only parton distribution in initial parton 
cascade is relevant for the large b t behavior. 

The above statement sounds strange because we claim that in the BFKL ladder we need to 
take into account the non-perturbative QCD correction only for two t-channel gluons attached 
to the larger dipole. Indeed, the large bt behavior stems from the exchange the lightest particles 
[1H] ( gluon in our case) somewhere in the middle of the 'ladder' (see Fig. EH)- Since in the 
middle of the ladder we have massless gluons the large bt behavior is determined by the diagram 
of Fig. EH It has a form 

N(k,k',y;b t ) = (7.119) 

-| J d 2 q dyi J d 2 (b t - b lt )N{k,q,y - yr,b t - b lt ) dy 2 J d 2 (b t - b lt )N(k,q,y 2 ;b 2t ) 

where we consider the bt is much larger than typical impact parameters in the integration over 
(b t — bit) and b 2t . Using Eq. (|fi.l()7j) we can take integrals over impact parameters in Eq. ()7.119j) . 
For the upper amplitude we have to consider S(b t ) = 1 in Eq. (|6.107|) since we do not introduce 
non-perturbative corrections, while for the lower amplitude we substitute our solution at large 
b t for which S(b t ) — ► e ~ 2m ^ b t. Finally, we have 

N(k,k',y;b t ) - b2 — j Q d yi F(y- yi ) ^ dy 2 y\. (7.120) 

Since F(y — y{) behaves as an exponent the integral chooses small y 2 . Therefore, actually the 
cell with large b t is the closest one to the largest dipole. This simple estimates show that our 
statement that the large bt behavior of the solution to the non-linear equation is determined 
by the initial condition is self-consistent . 

The great advantage of the semi-classical approach is that we clearly see all phases of high 
parton density QCD (see Fig. I2U|) . The critical line separates the phase with high parton 
densities. To the left of critical line we have a phase in which we have both saturation and 
the geometrical scaling behavior for the gluon density. This phase L.McLerran and his team 
suggested to call "Color Glass Condensate" (CGC) since the condensed system of color dipoles 
behaves as glass : a disordered system which evolves very slowly relative to a natural time 
scale (see more and better in a new review of this approach written by Iancu and Venogapalan 
in Ref. ^!)]).To the right of the critical line the parton system is still rather dense but it is 
far away from saturation while still has a geometrical scaling behavior [20] ■ This system be 
can call Color Dipole Liquid (CDL). The order parameter for CGC — > CDL transition is the 
saturation scale given by the critical line. The trajectory which is tangent to the critical line, 
separates the CDL phase from the phase with low density of color dipoles for which we can use 
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Figure 19: The diagram which de- 
scribes the large b t behavior of the 
dipole amplitude. Two gluons are 
close to the mass shell. 



the linear evolution equations: the BFKL or DGLAP ones. The natural name for this phase is 
Color Dipole Gas. 
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